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We report two novel determinations of \V u b\ with reduced model dependence, based on mea- 
surements of the mass distribution of the hadronic system in semileptonic B decays. Events 
are selected by fully reconstructing the decay of one B meson and identifying a charged lepton 
from the decay of the other B meson from Y(45) — > BB events. In one approach, we combine 
the inclusive B — > X u iv rate, integrated up to a maximum hadronic mass nix < 1.67GeV/c 2 , 
with a measurement of the inclusive B — > X s y photon energy spectrum. We obtain \V u b\ = 
(4.43 ± 0.38 sta t ± 0.25 syst ± 0.29 the o) x 10~ 3 . In another approach we measure the total B — > X u £f 
rate over the full phase space and find \V u b\ = (3.84 ± 0.70 s tat ± 0.30 sys t ± O.lOthoo) x 10 -3 . 

PACS numbers: 13.20.He, 12.15.Hh, 14.40.Nd 



The measurement of the element V U b of the Cabibbo- 
Kobayashi-Maskawa quark- mixing matrix plays a crit- 
ical role in testing the consistency of the Standard Model 
description of CP violation. The uncertainties in existing 
measurements 0, Q are dominantly due to uncertainties 
in the 6-quark mass mi and the modeling of the Fermi 
motion of the b quark inside the B meson In this 
paper, we present two techniques to extract \V u b\ from 
inclusive B — > X u lv || decays where these uncertain- 
ties are significantly reduced. Neither method has been 
previously implemented experimentally. 

Leibovich, Low, and Rothstein (LLR) have presented 
a prescription to extract | V u b | with reduced model depen- 
dence from either the lepton energy or the hadronic mass 
mx 0- A technique utilizing weight functions had been 
proposed previously by Neubert Q. The calculations 
of LLR are accurate up to corrections of order a 2 and 
(Atob/(^to;,)) 2 , where £ is the experimental maximum 
hadronic mass up to which the B — > X u iv decay rate is 
determined and A « Aqcd- This method combines the 
hadronic mass spectrum, integrated below £, with the 
high-energy end of the measured differential B — > X s j 
photon energy spectrum via the calculations of LLR. 

An alternative method 7] to reduce the model de- 
pendence is to measure the B — > X u £v rate over the 
entire mx spectrum. Since no extrapolation is nec- 
essary to obtain the full rate, systematic uncertainties 
from rrib and Fermi motion are much reduced. Perturba- 
tive corrections are known to order a 2 . We extract the 
B — > X u lv rate from the hadronic mass spectrum up to 
C = 2.5GeV/c 2 which corresponds to about 96% of the 
simulated hadronic mass spectrum. 

The measurements presented here are based on a 
sample of 88.9 million BB pairs collected near the 
T(4S) resonance by the BABAR detector || at the PEP- 
II asymmetric-energy e + e~ storage rings operating at 
SLAC. The analysis uses T(4S) — * BB events in which 
one of the B mesons decays hadronically and is fully re- 
constructed (B r ) and the other decays semileptonically 
(B s \). To reconstruct a large sample of B mesons, we fol- 



low the procedure described in Ref. [2j in which charged 
and neutral hadrons are combined with an exclusively 
reconstructed D meson to obtain combinations with an 
energy consistent with a B meson. While this approach 
results in a low overall event selection efficiency, it allows 
for the precise determination of the momentum, charge, 
and flavor of the B r candidates. 

We use Monte Carlo (MC) simulations of the BABAR 
detector based on GEANT4 to optimize selection crite- 
ria and to determine signal efficiencies and background 
distributions. Charmless semileptonic B — > X u lv de- 
cays are simulated as a combination of resonant three- 
body decays (X u = tt, p, oj, 77, r/) Il0|. and decays to 
non-resonant hadronic final states X u [llj for which the 
hadronization is performed by JET SET 7.4 The ef- 

fect of Fermi motion is implemented in the simulation 
using an exponential function with the parameters 
m b = 4.79GeV/c 2 _and Ai = -0.24GeV 2 /c 4 13]. The 
simulation of the B — > X c lv background uses a Heavy 
Quark Effective Theory parameterization of form factors 
for B-^ D*tD and models for B -> Dit£v, D*ir£v 
and B -> D£v, D**£v [Hj decays. 

Semileptonic B s i candidates are identified by the pres- 
ence of at least one electron or muon with momentum 
p} > lGeV/c in the B s \ rest frame. For charged B T 
candidates, we require the charge of the lepton to be 
consistent with a primary decay of a B s \. For neutral 
B r candidates, both charge-flavor combinations are re- 
tained and the average B°-B° mixing rate ,16| is used to 
determine the primary lepton yield. Electrons (muons) 
are identified [l7| (Ref. @), with a 92% (60-75%) aver- 
age efficiency and a hadron misidentification rate ranging 
between 0.05% and 0.1% (1-3%). 

The hadronic system X in the B — > Xlv decays is re- 
constructed from charged tracks and energy depositions 
in the calorimeter that are not associated with the B T 
candidate or the identified lepton. The neutrino four- 
momentum p v is estimated from the missing momentum 
four- vector p miss = p T (4S) - Pb, ~Px ~Pe, where all mo- 
menta are measured in the laboratory frame and px{4S) 



5 






• Data 

u 1 v 


i 


' ! ' 1 



In terms of measurable quantities, S1Z U (£) is 



_sl -si 



3 4 5 
m x [GeV/c 2 ] 



FIG. 1: The mx distributions (without combinatorial back- 
grounds) for B — -> X£i> candidates: a) data (points) and fit 
components after the minimum-x 2 fit, and b) data and sig- 
nal MC after subtraction of the B — ► X c lv and other back- 
grounds. The upper edge of the eighth bin is chosen to be 
at m x = 2.5 GeV/ c 2 . This fit result, with X 2 = W.2 for ff 
degrees of freedom, is used to extract the number of signal 
events below 2.5GeV/c 2 . 



refers to the T(4S') momentum. 

To select B — » X u lD candidates we require exactly one 
lcpton with pg > 1 GeV/c in the event, charge conserva- 
tion {Qx + Qe + QBr = 0), and a missing four- momentum 
consistent with a neutrino hypothesis, i.e., missing mass 
consistent with zero (—1.0 < TO 2 liss < 0.5 GeV 2 /c 4 ), 
IPmissI > 0.3GeV/c, and | cos 6» miss | < 0.95, where 9 miss 
is the polar angle of the missing momentum three- vector 
p miss . These criteria suppress the majority of B — > X c iD 
decays that contain additional neutrinos or an unde- 
tected K® meson. Additionally we reject events with 
charged or neutral kaons (reconstructed as K® — > 7r + 7r~ 
decays) in the decay products of the B s \. We sup- 
press B — > D*iu backgrounds by partial reconstruction 
of charged and neutral D* mesons via identification of 
charged and neutral slow pions. The reconstruction of 
the mass of the hadronic system is improved by a kine- 
matic fit that imposes four-momentum conservation, the 
equality of the masses of the two B mesons, and p 2 — 0. 
The resulting mx resolution is ~ 250MeV/c 2 on average. 

The extraction of | V^b| /J | from the selected events 
starts from the equation 6] 

V*\ f6a(l + ^)(Cf))» x 1 



IV, 



7r[/ (C)+/+(C)] 



where 51Z U (£) is the partial charmless semileptonic decay 
rate extracted from the number of B — > X u lD events up 
to a limit £ in the mx spectrum. H^ax accounts for 
interferences between electromagnetic penguin operator 
Oj with O2 and Os an d £7°^ is the effective Wilson 
coefficient. The terms io(C) an d ^+(C) are determined by 
multiplying the photon energy spectrum dT 1 /dE 1 inB— > 
X s j decays 0] with weight functions and integrating. 
The weights are zero below a minimum photon energy 
E™ in = m B /2 - C/4. 



(2) 



Here, N U {Q is the number of reconstructed B — > X u lv 
events with mx < d /(C) accounts for migration in 
and out of the region below £ due to finite mx reso- 
lution, B(B — > XiD) is the total inclusive semileptonic 
3 4 5 branching fraction, and £ V (Q is the efficiency for select- 
nix [GeV/c ] i n g B — > X u iD decays once aB-> XiD decay has been 
identified with a hadronic mass below £. N s \ is the num- 
ber of observed fully reconstructed B meson decays with 
a charged lepton with momentum above 1 GeV/ c, ef/ef 
corrects for the difference in the efficiency of the lep- 
ton momentum selection for B — > Xlv and B — > XJLv 
decays, and £reco/ e reco accounts for the difference in 
the efficiency of reconstructing a B r in events with a 
B — > XiD and B — ► X u iD decay. By measuring the ratio 
of B — » X u iD events to all semileptonic B decays many 
systematic uncertainties cancel out. 

We derive N u (£) from the mx distribution with a 
binned \ 2 n t to four components: data, B — ► X u lD 
signal MC, B — > X c iD background MC, and a small 
MC background from other sources (misidentified lep- 
tons, B — > XtD t , and charm decays), fixed relative to 
the B — > X c iD component. N U (Q is determined after 
the subtraction of the fitted background contributions. 
For all four contributions, the combinatorial background 
is determined, separately in each bin of the mx distri- 
bution, with unbinned maximum likelihood fits to dis- 
tributions of the beam energy-substituted mass Toes = 
y/s/4 — pj of the B r candidate, where y/s is the e + e~ 
center-of-mass energy. The toes fit uses an empirical de- 
scription of the combinatorial background shape [19| with 
a signal shape J20| peaking at the B meson mass. The 
combinatorial background varies from 5% (low mx bins) 
to 25% (high mx bins). The fitted mx distributions are 
shown in Fig. ^ before (a) and after (b) subtraction of 
backgrounds. The mx bins are 300MeV/c 2 wide except 
that one bin is widened such that its upper edge is at £■ 
We extract V sl = (3.253 ± 0.024) x 10 4 from an un- 
binned maximum likelihood fit to the toes distribution 
of all events with p* g > 1 GeV/c. The efficiency correc- 
ts) tions ef/e l l = 0.82 ± 0.02 stat , as well as e u (C) and /(C) 
(see Tablefl) are derived from simulations, where we also 
find £re CO /e" eco in agreement with one, assigning a 3% 
uncertainty. 

We study three categories of systematic uncertainties 
in the determination of \V u b\- uncertainties in the signal 
extraction, the simulation of physics processes, and the 
theoretical description. The quoted uncertainties have 
been determined for a value of C = 1-67 GeV/ c 2 where 
the total uncertainty on \V u b\ is found to be minimal. 

Experimental uncertainties in the signal extraction 
arise from imperfect description of data by the detector 
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TABLE I: Quantities in Eq.[5]that depend on £ and their sta- 
tistical uncertainties. The LLR (full rate) technique is given 
in the first (second) column. 



c 


1.67GeV/c 2 


2.50 GeV/c 2 


/ 


f.010 ± 0.005 


0.998 ± 0.002 


N u 


120 ± 17 


135 ± 45 


Eu 


0.231 ± 0.005 


0.231 ± 0.004 


5TZ U x 10 3 


1.43 ± 0.21 


1.59 ± 0.53 





L LLR 


Full 


\ ! 




Rate 

















1 2 3 

£ [GeV/c 2 ] 



FIG. 2: \V u b\ as a function of £ with the LLR method (left) 
and for the determination with the full rate measurement 
(right). The error bars indicate the statistical uncertainty. 
They are correlated between the points and get larger for 
larger £ due to larger background from B — ► X c lv. The total 
shaded area illustrates the theoretical uncertainty; the inner 
light shaded (yellow) area indicates the perturbative share of 
the uncertainty. The arrow indicates £ = 1.67GeV/c 2 . 

simulation. We assign 0.5% (0.5%, 0.8%) for the particle 
identification of electrons (/i, K^), 0.7% for the recon- 
struction efficiency of charged particles, and 0.8% for the 
resolution and reconstruction efficiency of neutral parti- 
cles. An additional 0.9% uncertainty is due to imperfect 
simulation of interactions. By changing the func- 
tion describing the signal shape in tties to a Gaussian 
function and switching from an unbinned to a binned 
fit method we derive an uncertainty of 2.2%. An uncer- 
tainty of 0.8% is determined by letting the contribution 
from other sources (see above) to the mx spectrum float 
freely in the minimum-^ 2 fit. The uncertainties on the 
inclusive B — > X s j photon energy spectrum are propa- 
gated including the full correlation matrix between the 
individual bins. 

The second category of systematic uncertainties arises 
from imperfections in the composition and dynamics 
of decays in the simulation, both in signal and back- 
ground. The uncertainties in the branching fractions of 
B -► D(*<**HuX decays contribute 0.7%. The un- 
certainties in the form factors in B — > D*lv decays 01 
introduce a 0.3% uncertainty. Branching fractions of D- 
meson decay channels contribute 0.2%. The relative 
contribution of the non-resonant final states has been var- 
ied by 20% resulting in an uncertainty of 0.5%. The 



TABLE II: Summary of results and uncertainties on \V u b\ for 
both approaches. The LLR (full rate) technique is given in 
the first (second) column. 





1 f\7 
l.D / 


9 ^ 


IK J x 10 3 


4.43 


3.84 


B — > X u lv stat. 


7.7% 


18.2% 


experimental syst. 


3.3% 


3.6% 


background model 


1.0% 


3.8% 


signal model 


3.9% 


5.6% 


theoretical 


6.2% 


2.6% 


B^X S "/ (stat., syst.) 


3.5%, 2.0% 




\V c b\ (exp., theo.) 


1.0%, 1.7% 





branching fractions of the resonant final states have been 
varied by ±30% (tt, p), ±40% (w), and ±100% (r) and 
rf simultaneously) resulting in an uncertainties of 1.0%. 
An uncertainty of 0.7% due to imperfect description of 
hadronization is determined from the change observed 
when we saturate the spectrum with the non-resonant 
component alone. We derive a 1.3% uncertainty due to 
the imperfect modeling of the KK content in the X u 
system by varying the fraction of decays to ss-pairs by 
30% for the non-resonant contribution |2l| . Even though 
the extraction of \V u b\ does not explicitly depend on a 
model for Fermi motion, there is still a residual depen- 
dency via the simulation of signal events. By varying the 
Fermi motion parameters nib and Ai within their respec- 
tive uncertainties, taking correlations into account |l3| . 
we derive an uncertainty of 3.5%. 

We calculate theoretical uncertainties in the weighting 
technique by varying the input parameters and repeating 
the weighting procedure including the calculation of all 
variables: H^ ix , as, and Wilson-coefficients. We vary 
a between a(m(,) and a(myy) with a central value of 
1/130.3 and find an uncertainty of less than 1%. For 
perturbative effects, an uncertainty of 2.9% is derived by 
varying the renormalization scale /i between nib/ 2 and 
2nib- Non-perturbative effects are expected to be of the 
order (Am B /(Cm t )) 2 , where A = 500MeV/c 2 [12, re- 
sulting in an uncertainty of 5.4%. Theoretical uncer- 
tainties in the measurement via the full rate are taken 
from Ref. 23| to be 1.2% (QCD) and 2.2% (HQE). 
Table [H] provides a summary of the uncertainties for 
C = 1.67GeV/c 2 and for C = 2.5 GeV/c 2 . 

Finally, we present two different determinations of 
\V u b\- First, using the weighting technique with the pho- 
ton energy spectrum in B^X s -f decays from Ref. [13|, 
the hadronic mass spectrum up to a value of £ = 
1.67GeV/c 2 , we find \V ub \/\V ts \ = 0.107 ± 0.009 stat ± 
0.006 sys t ± 0.007thco- If we assume the CKM matrix is 
unitary then \V ts \ = \V cb \ x (1±C(1%)) and, taking \V cb \ 
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from Ref. [24|. we derive 

\V ub \ = (4.43 ± 0.38 ± 0.25 ± 0.29) x 10~ 3 , 

where the first error is the statistical uncertainty from 
B — > X u £v and from B — > X s j added in quadrature, 
the second (third) is systematic (theoretical). Second, 
we determine \V u b\ from a measurement of the full mx 
spectrum, i.e., up to a value of £ = 2.5 GeV/c 2 , and find 
\V ub \ = (3.84±0.70 st at±0.30 syst ±0.10 thco ) x 10" 3 usin 



the average B lifetime of t b = (1.604± 0.012) ps |l&|2i 

The weighting technique is expected to break down at 
low values of £, since only a small fraction of the phase 
space is used. Figure |5] illustrates the dependence of the 
result, and its statistical and theoretical uncertainties, 
on variations of £ and also compares it with the value 
of \V u b\ determined from the full rate. The weighting 
technique appears to be stable down to ( ~ 1.4GeV/c 2 . 
The current uncertainties on the B—fX s j photon energy 
spectrum limit the sensitivity with which the behavior at 
high £ can be probed. 

The above results are consistent with previous mea- 
surements [HQ but have substantially smaller uncertain- 
ties from nib and the modeling of Fermi motion. Both 
techniques are based on theoretical calculations that are 
distinct from other calculations normally employed to ex- 
tract \V u b\ and, thus, provide a complementary determi- 
nation of | V u b | ■ 
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